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Abstract 

We emphasize that the stabihzing symmetry for dark matter (DM) particles does not have to be the com- 
monly used parity {Z2) symmetry. We therefore examine the potential of the colliders to distinguish models 
with parity stabilized DM from models in which the DM is stabilized by other symmetries. We often take 
the latter to be a Z3 symmetry for illustration. We focus on signatures where a single particle, charged under 
the DM stabilization symmetry decays into the DM and Standard Model (SM) particles. Such a Za-charged 
"mother" particle can decay into one or two DM particles along with the same SM particles. This can be 
contrasted with the decay of a Z2-charged mother particle, where only one DM particle appears. Thus, if the 
intermediate particles in these decay chains are Oj(7-shell, then the reconstructed invariant mass of the SM 
particles exhibits two kinematic edges for the Z3 case but only one for the Z2 case. For the case of on-shell 
intermediate particles, distinguishing the two symmetries requires more than the kinematic edges. In this 
case, we note that certain decay chain "topologies" of the mother particle which are present for the Z3 case 
(but absent for the Z2 case) generate a "cusp" in the invariant mass distribution of the SM particles. We 
demonstrate that this cusp is generally invariant of the various spin configurations. We further apply these 
techniques within the context of explicit models. 



1 Introduction 



There is compelling evidence for the existence of dark matter (DM) in the universe |1J. These 
observations can be explained by the postulating of new stable particles. A consensus picture of 
the nature of such a particle is provided by a host of astrophysical, cosmological and direct detection 
experiments: A viable DM candidate must be electrically neutral and colorless, non-relativistic at 
redshifts of z ~ 3000 and generate the measured relic abundance of /i^ ^bm = 0.1131 it 0.0034 [2]. 
Additionally a Weakly Interacting Massive Particle (WIMP) is a very well-motivated paradigm 
[T]. Consider DM particles as relics which were once in thermal equilibrium with the rest of the 
universe. It is well-known that the measured relic abundance is correlated with the dark matter 
annihilation cross section [3j by 

,2o O.lpb-c 

h^UBM ^ -, ^ • (1) 

[av) 

The annihilation cross section of a pair of dark matter particles into a two particle final state goes 
as ^ 

where g denotes the couplings and M the masses of the particles in the dark sector. This cross- 
section is naturally of the right value for g ~ 0{1) and M ~ 100 GeV. Moreover, many extensions 
of the SM at the weak scale, most of which are invoked primarily as solutions to other problems of 
the SM (most notably the Planck- weak hierarchy problem), contain such stable WIMPs. Because 
of this possibility, it may be possible to detect DM directly via scattering off nuclei or indirectly 
via detection of its (SM) annihilation products [l]. 

Such a scenario also makes the idea of dark matter amenable for testing at the high-energy 
colliders. It is possible to produce only DM particles directly at colliders, but then we do not have 
any visible signal since the DM particles will simply escape these detectors without interacting. 
Instead we investigate events where the dark matter is produced (indirectly) along with visible SM 
particles from the decays of particles charged under both dark matter stabilization (but heavier 
than the DM) and the unbroken SM symmetries. The existence of such "mother" particles is a 
feature of almost all models of physics beyond the SM that contain stable WIMPs. 

To date, a tremendous amount of effort has been made to reconstruct such events at the upcom- 
ing Large Hadron Collider (LHC) in order to determine the masses of the DM, the mother particles 
and possibly intermediate particles in the decay chains. For example, see references [U [5l El [7] . 
Most of this work has been for parity (Z2) stabilized dark matter. This is because the most popular 
models, e.g. super symmetric (SUSY), little Higgs and extra-dimensional scenarios [8l [9l [TO l [TT t [T2] . 
all ensure the dark matter candidates remain stable by employing a Z2 stabilization symmetry. 
Importantly these models have served as a guide of expected signatures of dark matter at the LHC 



2 



In this paper we emphasize that any discrete or continuous global symmetry can be used to 
stabilize dark matter^ Furthermore, because all fundamental particles in nature are defined by 
how they transform under various symmetries, most of the popular (Z2) models actually consider 
only one type of DM candidate! It is therefore critical to determine experimentally, i.e., without 
theoretical bias, the nature of the symmetry that stabilizes dark matter. We embark on a program 
of study to distinguish models in which the DM is stabilized by a Z2 discrete symmetry from models 
in which the DM is stabilized with other symmetries. A beginning effort was made in reference [16] 
which focused on signatures with long-lived mother particles, i.e., which decay to DM and the SM 
particles outside of the detector. In this paper we study the complementary possibility of mother 
particles which decay to DM and the SM particles inside of the detector. 

Our main idea is that the final states and the "topology" of the decay of a mother particle are (in 
part) determined by the DM stabilization symmetry. Thus reconstructing the visible parts of these 
decay chains will allow us to differentiate a model of DM stabilized with a mm-Z^ symmetry from 
one where DM is stabilized with a Z2 symmetry. In this paper we begin to explore such signatures. 
Our conclusions seem generic for most stabilization symmetries that are not parity symmetries; 
however, for definiteness, we focus on the case of a Z3 symmetry. When illustrating the signatures 
we will generically refer to any model stabilized with Zi and Z3 stabilization symmetry simply as 
Z2 and Z3 models, respectively. 

More specifically to see differences between Z2 and Z3 models, we focus on the kinematic edges 
and shapes of invariant mass distributions of the SM particles resulting from the decay of a single 
mother particle charged under the SM and the DM stabilization symmetries. We note the possibility 
of one or two DM particles in each decay chain being allowed by the Z^ symmetry (along with SM 
particles which can, in general, be different in the two decay chains). Whereas, in Z2 models, decays 
of a mother particle in given SM final state cannot have two DM particles in the decay chain and 
hence typically has only one DM particle. Thus, 

• If all the intermediate particles in the two decay chains are off-shell and the SM particles in 
the two decay chains are the same, then we show that there are two Z3 kinematic edges in 
the invariant mass distribution of this SM final state at approximately Mmother ~ "i-dm and 
Mjnother — 2mDM- Modcls with Z2 stabilized dark matter have only one endpoint approxi- 
mately given by Mmother - "t-dm- 

In the case of on-shell intermediate particles, the decay of such a mother in a Z3 model can 
similarly result in double edges due to the presence of one or two DM in the final state. However, 
in this case the endpoint also depends on the masses of intermediate particles. Thus it is possible 
to obtain multiple edges even from decay of a single mother particle in a Z2 model due to different 
intermediate particles to the same final state. Hence, multiple edges are not a robust way to 

^ Gauge symmetries alone cannot be used to stabilize dark matter. See the discussion in reference [15] . 
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distinguish between Z3 and Z2 symmetries in the case of on-sheh intermediate particles. For the 
case of on-shell intermediate particles, we thus use shapes of invariant mass distributions instead 
of edges. In particular, 

• We find a unique decay chain topology with two SM particles separated by a DM particle 
(along with another DM at the end of the decay chain) which is generally present for Z3 
models but absent for the Z2 case. Based on pure kinematics/phase space, this topology 
leads to a "cusp" (i.e., derivative discontinuity) in the invariant mass distribution of the SM 
particles. 

Of course for a generic model, it is possible to have a "hybrid" scenario where elements from the 
on-shell and off-shell scenarios are present. 

An outline of the paper is as follows: In the next section, we begin with the case of off-shell 
intermediate particles in a decay chain. There we show how differing kinematic edges can be used to 
distinguish Z2 from Z3 models. In section [3] we move on to the case of intermediate particles being 
on-shell. There we show the existence of a "cusp" in Z3 models for certain topologies; further we 
show that this cuspy feature survives even when taking spin correlations into consideration. We then 
discuss a couple of explicit models - one based on warped extra dimensional framework [n\ [T8] and 
another using DM stabilization symmetry from spontaneous breaking [15]; see also reference |19] 
for another example of a Z3-model. Here DM is not stabilized by Z2 symmetries. In the second 
model, we show our signal is invariant under basic detector and background cuts. We next conclude 
and briefly enumerate how Z2 models can fake signals from Z3 models. We also mention future 
work to better reconstruct and distinguish Z3 from Z2 models, e.g., using the two such decay chains 
present in each full event. 

2 OfF-shell intermediate particles 

In this paper, we mostly study the decay of a single heavy particle, which is charged under the dark 
matter stabilization and SM symmetries, into SM and dark matter candidate(s) inside the detector. 
Henceforth, we denote such heavy particles by "mother" particles. In this section we assume that 
all intermediate particles (if any) in this decay chain are Ojff -shell. This off-shell scenario has been 
frequently studied by the ATLAS and CMS collaborations \13\ [H] for SUSY theories (which is an 
example of a Z2 model). 

We consider constructing the invariant mass distribution of the (visible) decay products. Unlike 
for the Z2 case, for Z^ models a mother particle A can decay into one or two DM particles along 
with (in general different) SM particles. We mostly assume, just for simplicity, that there exist 
two visible particles (a, b or c, d) in the final state as shown below (note however that the same 
argument is relevant to the general cases where more than two visible particles are emitted)§]: 

^See Figure [9] for appearance of these two types of decays, including the required interactions, in the context of 
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b 
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(3) 



DM 



Here (and henceforth) the "blob" denotes intermediate particles in the decay which are off-shell. 
Also, upper-case letters/red/solid lines denote particles charged under the DM symmetry (Z3 or 
Z2) and lower-case letters/black/solid lines denote SM (or "visible", as opposed to DM) particles, 
including, for example, a W boson. Such an unstable SM particle decays further into SM fermions, 
at least some of which are observed by the particle detector. 

In order to avoid any possible confusion, we would like to explain the above diagrams. Under 
the Z3 symmetry, a particle/field (p transforms as 



where q = (i.e., Za-neutral) or g = 1,2 (non-trivial Zs-charge). Suppose the lightest of the Z3- 
charged particles (labeled 0o) has charge q = 1 (similar argument goes through for charge q = 2 
for 00 )• Clearly, its anti-particle (0o) has (a different) charge q = —1 (which is equivalent to q = 2) 
and has same mass as (pQ. Then, solely based on Z3 symmetry considerations, all other (heavier) 
Zs-charged particles can decay into this lightest Zs-charged particle (in addition to Za-neutral 
particles, including SM ones). To be explicit, a heavier Za-charged particle with charge q = i can 
decay into either (single) (po or two (po^s (and Za-neutral particles). Taking the CP conjugate of 
the preceding statement, we see that a heavier Zs-charged particle with the other type of charge, 
namely q = 2, is allowed to decay into two (po's or single (pQ. Of course, (po cannot decay and thus is 
the (single) DM candidate in this theory. We denote this DM particle and its anti-particle by DM 
and DM, respectively, in above diagram (and henceforth), although we do not make this distinction 
in the text since DM and anti-DM particles are still degenerateo 

For simplicity, we assume that the SM (or visible) parts of the event can be completely re- 
constructed □ Considering the invariant masses niab and rricd, which are formed by the two SM 
particles a, b and c, d in each decay chain, one can easily derive the minimum and the maximum 

an explicit model. 

^Of course, which of the two particles is denoted anti-DM is a matter of convention. Also, as a corollary, the DM 
particle should be Dirac fermion or complex scalar in a Z3 model. 

^We explore the effects of basic background and detector cuts at the LHC for a simple model in section 14.1.21 
There we show the effects discussed in this section remain after cuts for the background. 
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kinematic endpoints of the distributions of niab and rricd which are given by |20j : 



m 



max 
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,max 



m 



cd 



nia + mb, 

Mraother - "^DM {l^eit process of Eq. (I9D^ , 

nic + md, 

Mmother " 2 TTiDM (^Right process of Eq. % 



(5) 
(6) 
(7) 
(8) 



Physically, the lower limit corresponds to the case when the two visible particles a, b (and similarly 
c, d) are at rest in their center-of-mass frame so that they move with the same velocity in any 
Lorentz frame. The upper limit corresponds to the case in which the DM particle(s) are at rest in 
the overall center-of-mass frame of the final state. Both maxima are independent of the masses of 
the virtual intermediate particles. The point is that the upper endpoints in the two distibutions 
are different. 

2.1 Double edge 

An expecially striking/interesing case is when the SM particles in the two decay chains are identical: 




DM 




(9) 



As we show below, it is possible to obtain a double edge in the distribution of this SM final state. 
We begin with presenting a basic idea of this phenomenon, before going on to more details. 

2.1.1 Basic Idea 

Taking into account the fact that the visible particles of both decays are the same and assuming that 
both subprocesses are allowed, the experimental distribution (l/E) dV /dniab will contain events of 
both processes. In such a combined distribution, clearly, the endpoint of Eq. ([5]) - denoted now 
by rn'J^^^ - will become an edge in the middle of the distribution, which along with the overall 
kinematic endpoint given by Eq. ([6]), will give rise to a double edge signal. Assuming the two edges 
are visible, it is interesting that we can determine both the DM and mother particle masses by 
simply inverting Eqs. ([6]) and ([8]): 



-^mother 
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Figure 1: Invariant mass distribution (l/F) dT/dmab for the processes of Eq. ([9]). The masses of 
the mother particle A and of the DM particles are niA = 800 GeV and ttidm = 300 GeV and the 
SM particles a and b are assumed to be massless. The solid and dashed curves on the left panel 
represent the distributions for the 3-body decay and the 4-body decay, respectively. On the right 
panel, blue/dashed (highest peaked), red/solid, and green/dot-dashed (lowest peaked) curves show 
the combined distributions with branching ratios of 3-body to 4-body given by 1:3, 1:1, and 3:1, 
respectively. 



In particular, the distance between the two edges is identified as the DM mass. 

In contrast to the cases just considered, in Z2 scenarios only one or three DM particles (i.e., 
not two) are allowed in a single decay chain due to .^2-charge conservation (unless the process is 
triggered with an uncharged mother particle [7J). Independently of phase-space considerations, we 
note that in Z2 models the decay chain with three DM particles should be highly suppressed with 
respect to the one DM case. The reason for such an expectation is that a decay with three DM 
in the final state requires a vertex with four (in general different) Zo-charged particles which is 
typically absent, at least at the renormalizable level in most modelsjj Therefore with only one 
possible decay process (in terms of the number of DM in the final state) we can only observe a 
single kinematic endpoint in the invariant mass distributions in a Z2 model. 



2.1.2 Details 

Of course the visibility of such a signal depends on the shapes of the distributions of each subprocess 
as well as their relative decay branching fractions. The solid curve and the dashed plot in the left 
panel of Figure [T] illustrate the generic shape of the distributions for the two processes of Eq. ([9]) 
based only on pure kinematics, i.e., no effects of matrix element and spin-correlations. (Such 



^Compare this situation to the Z3 case, where appearance of two DM in a decay chain comes from vertex with 
three Zs-charged particles which is more hkely to be present, especially at the renormalizable level. 
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Combined decays: Off— shell 
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Figure 2: Same as the right panel of Figured] but using a smaller DM mass, ttidm = 50 GeV. The 
edge in the middle of the distribution is no longer apparent. 

effects might be important and we will return to this issue in the context of specific models to show 
that multiple edges can still "survive" after taking these effects into consideration.) Because of 
the phase-space structure of the processes one realizes that the distribution in the case of 3-body 
decays is more "bent" towards the right (i.e., larger values of invariant mass) whereas for the 4- 
body decays the peak of the distribution leans more towards the left (i.e., smaller values of invariant 
mass). Because of this feature, the combination of the two distributions can give rise to two visible 
edges (as long as the relative branchings of the two decays are of comparable size). This is shown 
in the right panel of Figure [1] in which we show the combined invariant mass distribution of the 
two visible SM particles, for three different relative branching fractions of the two subprocesses. 
Based on the location of the edges in right panel of Figured] and Eqs. (|10p and (|lip . the mass of 
DM particle must be about 300 GeV and the mass of the mother particle must be about 800 GeV, 
which are of course the masses used in the example. 

Whether or not the double-edge signal is clear (and hence we can determine the DM and mother 
masses) also depends on the DM mass which must be relatively sizable compared to the mass of 
the mother particle. For example, if we take a DM mass of 50 GeV instead of 300 GeV that we 
assumed above, with the mother mass fixed at 800 GeV, we observe from Figure [2] that the plotted 
distribution does not provide a good measurement of Mmother and mDM- 

Let us return to the issue of the relative branching fraction for each subprocess. The decay 
into two DM particles should be generically phase-space suppressed relative to the decay into just 
one DM particle. So, based on pure phase-space suppression, the branching ratio of the decay into 
two DM might be much smaller than the decay into one DM (unlike what is chosen in the figures 
above). Hence, it might be difficult to observe a double-edge signal. However, in specific models 
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this suppression could be compensated by larger effective couplings so that the two decays have 
comparable branching ratio, and therefore, the double-edge is visible as in Figure [TJ 

In fact, another possibility is that the two decay chains for the Z3 case, i.e., with one and two 
DM particles, do not have identical SM final states, but there is some overlap between the two SM 
final states. For example, 




If we assume that particle c is (at least approximately) massless, then the maximum kinematic 
endpoint of ruab in the first of the above-given two reactions is still M^other — ^dm — ~ 
Mmother — mDM- In this situation both the reactions have 4-body final states and hence could be 
easily have comparable rates, at least based on phase-space (c./. Earlier we had 3-body vs 4-body 
by requiring the same two-body SM final state for the two reactions found in ([9])). On the other 
hand, although the two rates are now comparable, it might actually be harder to observe a double 
edge because the shape of the two individual distributions are both peaked towards the left (i.e., 
smaller values of invariant mass) and even if they have different end-points, the combined distribu- 
tion might not show as clearly a double edge as the earlier case where the two shapes are apparently 
distinct. 

2.2 Different Edges in Pair Production 

Finally, what if there are no common SM particles between the final states of the decay chains with 
one and two DM particles so that we do not obtain a double edge? In this case, one can consider an- 
other analysis, by making the further assumption that the same mother particle A is pair-produced 
in each event, and that the decay products of each A are now distinct and very light or massless, i.e., 
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Here we have chosen three SM particles (a, b, c) in the decay chain with one DM just so that 
both decay chains involve a 4-body final state. In this situation one can restrict to events with all 
five SM particles (a,...e) particles in the final stat^, but use both sides of the event, i.e., obtain 
the full invariant mass distribution of the visible particles of each (distinct) side. In the inter- 
pretation of these results in the context of a model, the difference between the endpoints of 
each separate distribution will give the dark matter mass, and like before, the mass of the mother 
particle A can be found using a combination of the two end-points, i.e., m dm = i^^^c^ ~ "^^e^ 
and M^other = 2m-- - m . 

3 On-shell Intermediate Particles 

In this section, we consider the case where the mother particle decays into SM and DM via in- 
termediate particles which are all on-shell. Again, like in section [2] all particles are assumed to 
decay inside the detector. In this case, the endpoints of invariant mass distributions will depend 
on the masses of these intermediate states as well as the masses of the mother and the final state 
particles. Both in the Z2 and Z3 cases there will be more possibilities for the upper endpoints 
because of the possibilities of "Multiple topologies" and "Different Intermediate Particles" (to be 
explained below) for the same visible final state. Since even for the Z2 case it is possible to obtain 
multiple edges, finding multiple edges is not anymore a robust discriminator between Zi and Z3 
unlike the off-shell decay case. We then discuss a topology of the decay chain which does allow us 
to distinguish between the two models. 

3.1 Additional sources of Multiple Edges 

Here we discuss how it is possible to obtain multiple edges even if we do not combine decays of 
mother particle into one and two DM particles. 

3.1.1 Multiple Topologies 

For Z3 models we can expect multiple endpoints from the decays of the same mother particle into 
a given SM final state by combining the two decay chains with one DM and two DM particles, 
respectively, just as in the case of the decays with off-shell intermediate particles. However, this is 
not the only way of obtaining multiple endpoints, i.e., such a combination of decay chains with one 
and two DM is not essential. The reason is that there are multiple possible topologies even with 
the completely identical final state if it contains two DM, due to the various possibilities for the 
locations of two DM particles relative to the other SM particles in a decay chain. For example, for 
the case of a 4-body decay process (i.e., two SM and two DM particles) there will be three different 
possibilities: 

®If we include other events which have a, 6, c or d, e on both sides, we still get the different edges that we discuss 
below, but as we will mention later, such events will not allow us to get rid of "faking" Z2 models. 



10 



DM a (^4^ 

D C B ~DM 

DM c a 

(15) 

D C B ~DM 

DM 

(16) 

D C B DM 

Note that (as above) decay cascades involve a "charged-charged-charged" (under Z3 symmetry) 
vertex (in addition to "charged-charged-neutral" vertices) in order to contain two DM particles in 
the final state. 

Assuming that the visible particles are massless, rua = rric = 0, the upper endpoints for each 
topology are given by (See Appendix A) for details.): 

= , , i'"h'"""::lr"'°'^ ^ > (!o.E^.m) 

(m--)2 = ("^C-"^I)("^|-^dm) I'forEq. (HSI)) (18) 

V / 

(m--)2 = ^^\-<i){fc-^\) fforEq. (USD) (19) 

where A is the well-known kinematic triangular function given in the form of 

A(x, y, z) = + + - 2xy - 2yz - 2zx. (20) 

The main point is that kinematic endpoints are functions of the masses of the mother, the DM and 
the intermediate particles, and moreover, this dependence changes according to different topologies. 
Thus, even if the intermediate particles involved in these decays of a given mother particle are the 
same, one will still obtain multiple endpointslll Finally, if we combine decay chains with one and 
two DM in the final state (even if the latter has just one topology), the difference between the two 
endpoints will not lead to a direct measurement of the DM mass like in the off-shell decay case 
because again, the mass of intermediate particles is one of the main ingredients to determine the 
endpoints. 

In Z2 models the decay topologies must have a single DM particle and that too at the end of the 
decay chain because the vertices in the decay cascade are of the form "odd-odd-even" (under the 
Z2 symmetry) Ifl Nevertheless, there can still be different topologies because of different ordering of 



^Of course, the different possible decay topologies can, in general, have different intermediate states. 
*Note that a similar argument applies to decay chains in Z3 models with only one DM in the final state. 
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the visible states. For example: 



c b a 




DC B DM 

b c a 




Obviously, the endpoints for a given invariant mass distribution, say rrica, will be different for each 
of these two topologies, and actually they can be obtained from Eqs. pT|) and ([T8|) by just replacing 
rriDM in the denominator of Eqs. (I17p by mf, and leaving Eqs. (llSp unchanged (and where and 
rric are still assumed to vanish). 

3.1.2 Different Intermediate Particles for Same Final State 

In addition, even if the topology and the order of visible particles are the same, there is the 
possibility of multiple paths for the same mother particle to decay into the same (SM and DM) final 
state by involving different intermediate particles. We will obtain multiple endpoints in this case 
because of the dependence of the endpoints on the masses of intermediate particles (as mentioned 
above). This argument is valid for both the Z2 and Z3 models (and one or two DM for the latter 
case): for a final state with two SM and one DM, we can have 



b a b a 




C B DM C B' DM 



For example, in SUSY, the decay chain X2 ^^^ Xi can proceed via intermediate right- or left- 
handed slepton. Since the masses of intermediate right- and left-handed sleptons are in general 
different, multiple endpoints are expected. 

3.2 Cusp Topology 

So far, we have learned that for on-shell intermediate particle cases the multiple edge signal is not 
a good criterion to distinguish Z3 from Z2. Instead, we focus on shapes of these distributions. 
Consider the topology which can be present in Z3 models (but absent in the Z2 case) with two 
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visible SM particles separated by a DM particlcl, i.e., 



DM 



D 



C 



B 



(24) 



A 



We assume massless SM particles (i.e., rua = rric = 0) and the mass hierarchy mo > mc > 

niB > ruA- Also, we neglect spin-correlation effects in this section. We sketch the derivation of the 

distribution of the ac invariant mass here and refer the reader to the Appendix A for details. The 

1 dV 

differential distribution — — — ^ that we want to study can be obtained for this "new" topology 

1 d'^T 

easily by noting that the differential distribution —— — — must be flat, where the variables are 



defined as follows 



r dudv 



1 — cos 9 



(C) 
cDM 



and 



1 — cos 6. 



{B) 
ca 



(25) 



with Oca being the angle between c and a in the rest frame of B, and 0).iyyi being the angle between 



c and DM in the rest frame of C [21]. In addition, we have Q < u,v < 1. Thus, we can write 

1 d'^T 



One further finds that 



r dudv 



e{i - u)e{u)e{i - v)e{v) 



(26) 



mla = mZ'^il - au)v, 



(27) 



where m™"^ is given in Eq. ()17p with raoM in the numerator replaced by m^, and so we can 
make a change of variables from the differential distribution of Eq. (|26p and obtain the distribution 

1 a^r 

— 7—— — which can then be integrated over u to finally obtain the distribution with respect to 
r dijkmnj.^ 
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r drric 
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< mc„ < \/r 
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mf. 



a m^^"^ < mr 



a m,. 



< m" 



(28) 



where m^^ is given in Eq. ()17p and 



a 



2Ai/2( 



m 



DM + Ai/2(?7i^, m^r^,m 



2 V 

dm) 



(29) 



®Note that in general D might come from the decay of another ^s-charged particle and similarly, at the end of 
the decay, A might not be the DM, that is, it could itself decay further into DM particles and other visible states as 
long as Za-charge conservation is respected. The to the left of D and to the right of A signify this possibility. 
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Note that 



1 dT 
Tdm 



2m 



1 dr 

r 9m2 
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4— body decay: Inv. mass distribution 4-body decay: Inv. mass squared distribution 




«iea (GeV) (GeV^) 



Figure 3: The panel on the left shows the distribution in rrica while the right hand panel shows 
the distribution in m^^ from the decay chain of Eq. (I24p . The masses of mother particle, two 
intermediate particles, and DM particles are 800 GeV, 700 GeV, 400 GeV, and 200 GeV, respectively 
and the SM particles are assumed massless. A "cusp" due to the topology of Eq. is clear in 
both distributions. 



From these results we can easily see that the new topology introduces two different regions in the 
nica distribution with a "cusp" at the boundary connecting both regions, located at \/l — a rn^^. 
Figure [3] shows the same distribution in both panels, but with respect to rrica on the left panel and 
with respect to m^^ on the right panel. As we will argue later, the second option seems better 
suited once spin correlations are taken into account, but in both plots, one observes that the cuspy 
feature is quite clear. 

3.2.1 Two Visible Particles 

Consider first the simple case of only two visible particles in a decay chain. In the Z3 reaction 
of Eq. (j24p . D is then the mother particle and A is DM. Clearly, we would find a cusp in the 
invariant mass distribution of the two visible particles in the Z3 model, but not for the Z2 model 
since the two visible particles must always be adjacent to each other in the latter case0 Thus, the 
presence/absence of cusp could be used to distinguish Z3 and from Z2 models. 

3.2.2 Generalization to More than Tw^o SM Particles in Decay Chain 

Of course, in general in both Z2 and Z3 models there will be more than two visible particles with 
possibly some of them being identical, and this will undoubtedly complicate the analysis. For 
example, in the reaction of Eq. ()24p . a, c, or both can be produced at some other place of the same 

^^Note that we are considering decay of a Z-i or Z2-charged mother. A Z2-uncharged, i.e., even, mother is allowed 
to decay into two DM and can give a cusp in the invariant mass of two visible particles from such a decay [?]■ 
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decay chain in addition to the locations shown there, e.g., 




D C B A DM 

Here a', which is an identical particle to a, is assumed to come from the immediate left of D, and 
c', which is an identical particle to c, is assumed to come from the immediate right of A. Note 
that there is no DM between a' and c (unlike between a and c) in first reaction above (similarly 
between a and c' vs between a and c in second reaction above), and that a and a', and c and c' 
are identical. Therefore, in both these examples, it is clear that we will obtain a more complicated 
distribution in mac than the one studied previously. 

Nevertheless, the method described previously to disentangle the Z2 from the Z3 cases (when 
having two visible particles), can still be generalized to the situation of many visible particles in a 
decay chain. For example, let us consider the case of three visible SM particles in the final state for 
both Z3 and Z2 models. We will obtain a cusp even in the Z2 case when considering the invariant 
mass of two not "next-door neighbor" visible particles such as in mac for the decay process in Eq. 
(pT|) . The reason is that, even though the precise topology of Eq. (p^ is absent in a Z2 model, a 
similar one is generated by the presence of a SM particle (i.e., b) in-between two other SM particles 
(i.e., a and c) as in Eq. ()2ip . Thus the analysis performed earlier for Eq. ()24p applies in this case, 
but with the DM mass set to zero (assuming SM particle b is massless). 

However, this type of degeneracy between Z2 and Z3 can be resolved by considering all of the 
three possible two- (visible) particle invariant mass distributions. In the Z3 case with two DM 
particles in the final state, two of these three invariant mass distributions will have cuspy features 
whereas only one such invariant mass distribution will have a cusp in the Z2 case. The reason is 
again that in the Z3 case, since one more particle is added to the decay products compared to the 
Z2 case (i.e., we have two invisible and three visible particles), there will be final state particles 
(visible or not) in-between the two visible particles for two of the three pairings. This feature 
remains true for more visible particles, i.e., in general we will obtain more cusps in the invariant 
mass distributions in a Z3 model than in a Z2 model. 

3.3 Spin Correlations 

Once spin correlations are involved, the derivative discontinuity (cusp) might appear unclear. Nev- 
ertheless, it may still be possible to distinguish a Z3 model from a Z2 model by employing the 
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Inv. mass distribution of a and c 




550- 



»i\a (GeV^) 

Figure 4: Invariant mass distribution of particles a and c, from the decay chain shown in Eq. (j24p . 
including spin correlations, and such that the intermediate particle C has spin 1 and the interme- 
diate particle B has spin 1/2, and the couplings are chiral. The "cusp" in this distribution appears 
more defined than in Figure [3] where spin correlations were not considered. 

fitting method which we will show in the rest of this section. The basic idea is that the distribution 
dT /dm^^ of three-body decays in Z2 (i.e., one DM particle and two visible particles) can (almost) 
always be fitted into a quadratic function in m^^, whereas the distribution of the new topology of 
cannot not be fitted into a single quadratic function, that is, two different functions are required 
for fitting each of the two sub-regions of the distribution. Let us see how this works for a Z2 model 
(i.e., one DM and two visible particles) and a Z3 model (i.e., two DM and two visible particles) in 
turn. 

3.3.1 Z2 case: 1 DM + 2 Visible 

We can again make use of the same angular variables considered earlier for the case of this 3- 
body decay cascade, shown for example in Eq. (j23p . According to the references |21j and [22], the 
normalized distribution including spin-correlations is given by 

^^^=e{t)9{l-t)f{t) (32) 

where again we have defined the variable t as 

1 - cos e[^^ , , 

t = (33) 

Here f{t) is a function of t and 6^^ is the angle between particles a and b of Eq. (p3|) in the 
rest frame of particle B. One then notes that mf^j = {m^^)'^ t which basically means that the 



16 



distribution with respect to the invariant mass m^^ (which is of our interest) is essentially the same 
as the one with respect to t above. This means that the distribution in will have the functional 
form /. According to the reference j23j, such spin correlation functions are just polynomials of 
cos (i.e., (1 — 2t)). Moreover, if we restrict our consideration to particles of spin-1 at most, the 
maximum order in t of the polynomial is two, which means that the most general form of / will be 

/(t) = Ci + C2t + C3t^ (34) 

In turn, the invariant mass distribution we are interested in must therefore take the form (in the 
region between the endpoints enforced by the ^-functions) 

^ ^ c'l + cWL + C3"^L- (35) 



Tdml 



With the experimental data we can construct the invariant mass distribution, and we will be able 

to determine the three constants c'l, c'2, and C3 by fitting into a parabola in the m^^ variable. In 

other words, for any 3-body decay chain, with or without spin-correlation, it is always possible to 

fit the invariant mass distribution ^ into a curve quadratic in m? . 

r dmf^ 



3.3.2 Z3 case: 2 DM + 2 Visible 

We now consider the new topology of Eq. (|24p including the possibility of spin correlations. As in 
Section 3.1, we use the same angular variables u and v. However, the normalized distribution with 
spin correlations become a little more complicated than before 

1 d^r 

r dudv 

where again 



e{v)9{l - u)g{u)e{v)e{l - v)h{v) (36) 



_ 1 - cos 6*^^^ _ 1 - cos 9^J^^ 
u = , V = . (37) 

Like in the previous section, g{u) and h{v) are spin-correlation functions (cf. g = h = 1 without 
spin correlation discussed earlier) and again the invariant mass squared is given by 

ml = {m^n\l-au)v. (38) 



where a is the same kinematical constant defined in Eq. (j29p . As in the analysis without spin 
correlations, the two types of ^-functions will split the entire region into two sub-regions, with a 
cusp at the separation point, whose location is independent of the spin correlation effects (since it 
depends on purely kinematical constants a and m™"^). But unlike the scalar case (i.e., with no 
spin correlations), we have now two functions g{u) and h{v) which can change the shape of the 
distribution and in principle affect the derivative discontinuity (the cusp). 
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In detail, by the chain rule the previous normalized distribution can be modified and partially 
integrated to obtain 

fd^~Jo - au) ^^""^ [{mZ^-y{l-au)) ^ ^ 

where 



Wmax = Max 



(40) 



The two possible choices in the definition of Umax above arise when integrating ^ dm^^du ^^^^ 
respect to u due, in turn, to the integration limits enforced by the 6 functions. This leads to two 
different regions for the differential distribution such that in the first sub-region, we have < rrica < 



\/l — a and Umax = while for the second region, we have \/l — a < rrica < "T-m^^ ^.nd 

limax = X ( 1 ~ T-^sSw ) [ZI]- So far, most of the steps are similar to the case of no spin correlations 
except for the presence of the factors of spin correlation functions, g and h. 

It would seem that we need to know the "precise form of g and h in order to proceed further, i.e., 
in order to perform the integration in Eq. (|39p . However, for the purpose of determining whether 
or not there is a cusp, we will show that it is sufficient to know the fact that those spin-correlation 
functions must be second order polynomials in their argument as mentioned in the analysis of the 
Z2 case. Using this fact we can write down the above integrand as 

^ -g{u)h[-^—\ = . .. t^+ ,, ^ ,n (M + M^) + -^(fe4 + M + M^) 



1 — au \1 — au J (1 — au)^ (1 — au)^ 1 — au 

+ (67 + 68*) + ^9(1 - an), (41) 

where we have introduced the same variable t = rn^^/ {m™^^)'^ used for the 3-body decays and where 
the kinematical constants 6j will depend on the specific nature of the couplings and particles in the 
decay chain (i.e., they must be calculated on a case by case basis). The terms of the integrand are 
organized as a power series in (1 — au) ~ instead of in u - because of the simplicity of the former 
form. Integrating then gives 



1 dT 
f~dt 



b[ + 6^t + 6^5*2 for < t < y/T^ 

(42) 

b'( + h'^t + hlt^ + {b'i + h'lt + hlt^) log t for ^/T^ < t < 1 



where again, the kinematical constants h[ and 6-' are specific to each situation. Thus, even with 
spin correlations, the functional dependence on t (oc m^^) is quite simple; however, the crucial point 
is that it is different for each sub-region of the distribution. In particular this simple dependence 
in the distribution of (and not rrica) suggests that it may be more appropriate to consider 
the distribution of rn^^ instead of the distribution of ruca- In Figured] we show the m^^ invariant 
mass distribution for the decay chain of Eq. (p^ . but in the special case where particle C has 
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spin 1 and the intermediate particle i? is a fermion, and some of the couphngs are chiral. We used 
Madgraph [21] to generate events taking the particles a and c to be massless and taking tudm = 100 
GeV. One can compare the shape of this distribution with the one from the right panel of Figure 
|3]and see that in this case, including the spin correlation makes the cusp even more apparent. 

One of the main differences between the two subregions is that the first one has no logarithmic 
dependence in t while the second (in general) does have it. Of course, from Eq. (j41|) . we see that 
this logarithmic term could be suppressed for the case 64 = 65 = 65 ~ 0. However, even in this 
special case we would still have to employ different sets of coefficients in the two sub-regions as 
follows. The functional forms in both the regions are now quadratic in t, i.e., 



1 dV 
f~dt 



67 + f(2-a) + (^ + 68) t + 



bi{2-a) 
2(l-a)2 



1-a 



b2 + h + 



(62 - b^ + bj-bs)t- (63 + ^8 + 



2 . 



< t < y/l^^ 



(43) 



a < i < 1 



Considering just the constant terms, we see that it is possible to obtain identical functions in the 
two regions only if a = 1 and 61 = 62 = 0. However, using Eq. (j20p and Eq. (j29p . it can be shown 
that a is always (strictly) less than 1. In other words, it is highly unlikely that the distribution in 
each sub-region can be fitted successfully to the same polynomial of order two in t; the cusp will 
thus survive even in this case. 

In figure [5] we show the distribution (again obtained with Madgraph [24]) for the same decay 
chain as in Figure [H but where the chiral structure of some couplings has been modified from 
before. We see that the cusp feature is now less apparent, but one also sees that a full fit to a 
polynomial of order two (left panel) is not as good as a multiple-region fit (right panel), where 
the first part of the distribution is fitted to a polynomial of order two (see first line of Eq. ()42|) ). 
and the right side of the distribution is fitted to the functional form (with a logarithm) given in 
the second line of Eq. 



4 Example Models 

4.1 Stabilization Symmetries from Spontaneous Symmetry Breaking 

The most popular models of physics beyond the SM focus on solving the weak-Planck hierarchy 
problem by adding new particles at the weak scale. Some of these new particles can be stable as a 
consequence of a discrete (often a parity) symmetry that is a part of (or imposed on) the theory. 
Thus these particles can have the correct thermal relic abundance to constitute dark matter, i.e., 
dark matter is then a "spin-off" of solving the hierarchy problem. It may be possible, however, 
that the question of the origin of dark matter is not rooted in first solving the hierarchy problem. 
In this case, it is the thermal relic density which "fixes" the mass of the dark matter particles to 
the weak scale. As well, it is also known that the dark matter and baryon densities today are close 
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Inv. mass distribution of a and c Inv. mass distribution of a and c 




(GeV^) m\ (GeV^) 

Figure 5: Invariant mass distribution of particles a and c, as in Figure HI but with different chiral 
couplings. The cusp position is less apparent in this case but one can see (left panel) that a fit to a 
polynomial of second order as shown in Eq . ([35|) is not very good (that is, the Z2 interpretation). 
On the right panel we show the same distribution, with a different fitting function for the left side 
of the distribution and the right side (see Eq. (jl2]) ). consistent with the existence of a cusp, i.e., 
the Z3 interpretation. 

in size 

r^DM ~ 4.70baryon (44) 

which provides a hint to a possible common origin which may have a solution at the weak scale. 

With this background, in reference |15j . a model was introduced in which a SU{N) gauge group 
is spontaneously broken to the Z^r center. There the goal was to, in part, determine whether a 
"copy" of weak interactions could generate a viable dark matter candidate. As is well known, the 
SM electroweak gauge group is spontaneously broken to ?7(l)em which makes the lightest electrically 
charged particle (i.e., the electron stable). By analogy the SU{N) gauge group in this model is 
broken to the Zi\j center that stabilizes dark matter. The dark matter candidates which transform, 
e.g., as a fundamental under the SU{N) are stabilized by the Zj\f. A unique, testable feature of 
these models is the existence of new gauge bosons that are neutral under the SM symmetries but 
do not couple to SM particles as Z bosons. These gauge bosons transform as adjoints under the 
SU{N) and therefore are invariant under the Zjv center. 

The fermionic dark matter candidate mentioned above (i.e., transforming as a fundamental 
under SU{N)) would get a mass m ~ Xv^ew Here A is a generic 0(1) Yukawa coupling. For this 
mass to be of order the weak scale as required for the thermal relic abundance [15], the vacuum 
expectation value (vev) that breaks the SU {N) must be of order the SM Higgs vev. Thus, the 
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SU{N) gauge bosons also have weak scale masses. As a bonus (or "double-duty") of the new 
SU{N) vev being similar to the Higgs vev, at finite temperature these gauge bosons have sphaleron 
solutions and nucleate additional bubbles (in analogy with the weak interactions) that may be 
relevant for electroweak baryogenesis. 

Let us make a connection to our study of distinguishing Z3 from Z2 using double edges described 
in section [2] and new topologies described in section O There, we chose symmetry mainly 
for illustration; in any case, the key to our signal is the existence of "charged-charged-charged" 
couplings under the dark matter stabilization symmetry in the decay chain. Parity stabilized 
models have only "odd-odd-even" couplings. Any model, not necessarily a Z3, that features such 
coupling has the potential to generate the signals we discussed earlier. 

To this end, we take a "toy" limit of the model discussed in reference |15] . Namely, we consider 
a scenario where the new gauge bosons are long-lived and register as missing energy (pr) in the 
detectors so that they behave effectively as dark matter particles, i.e., as "charged" (even though 
they were uncharged "to begin with"). This assumption will then "convert" the "dark" SU{N) 
gauge coupling of a SU{N) fundamental fermion into an effective "charged-charged-charged" cou- 
pling which will result in the double kinematic edges as well as the new topologies discussed earlier. 
The result will also be to generate a "hybrid" of the on- and off-shell scenarios presented above. 

Here, we simply want to make an estimate of the robustness of the signal described in sections [2] 
and [3] in the presence of basic detector and background cuts. So, the above toy limit of model 
in reference [15j will suffice for such a study of exploring the effects of the "charged-charged- 
charged" coupling in a more realistic situation than considered in earlier sections o As a first 
step, following |15) . we summarize the effective lagrangian for our model. Later we will discuss a 
simple production mechanism and discuss cuts consistent with the ATLAS and CMS collaborations. 
Results follow afterwards. 

For simplicity we chose = 2 to make our analysis. To break the SU{2)£) — t- Z2 we require 
two new additional Higgses in the "dark" sector which transform as an adjoint 

(45) 




The higgs generate the following vevs 



(46) 



which break the SU{2)u to the center. A general scalar potential does not naively generate the 
required breaking. To get the correct vacuum alignment, we require a scalar potential in the "dark" 



^^Alternatively, extensive model building along the lines of reference [TS], which is not the focus of this paper, can 
provide a "genuine", i.e., without assuming long-lived gauge bosons, "charged-charged-charged" coupling. 
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sector to minimizes cj) ■ rj. SU (2) ^ scalar potential is 

2 



V 



= Xii^(P'^ + X7(j)-r]-vlj +A2(r/2 + A8</>-7?-wi ) + A3 ( + r/^ 



(47) 



which generates three new heavy gauge bosons as well as three additional Higgses. In addition, we 
add anomaly free scalar and fermions with the quantum numbers listed in Table 1. Constructing a 



Particle 


SU{3)c 


SU{2)l 


SU{2)d 


U{1)y 


Q 


3 


1 


2 


1/3 


sQ 


3 


1 


2 


1/3 


L 


1 


2 


2 


-1/2 


sL 


1 


2 


2 


-1/2 


X 


1 


1 


2 





SX 


1 


1 


2 





y. 


1 


1 


3 







1 


1 


3 





V 


1 


1 


3 






Table 1: An effective, anomaly free particle spectrum that fills out a (5, 2)-|-(5, 2). The "s" prefactor 
denotes a scalar particle. Here are the SU{2)d gauge bosons. We assume the mass of the Q and 
SX is heavy and integrated out. The rest of the spectrum mediates the decay chain in equation! 



supersymmetry UV completion to this effective lagrangian is straightforward. Although the details 
is beyond the scope of this paper, note a simple way to do so would be to augment minimum su- 
persymmetric standard model with chiral superfields with the charges in Table ICT SUSY breaking 
terms would then need to be constructed to lift the appropriate particles which will be integrated 
out to generate the effective theory. 

4.1.1 Production Rates at the LHC 

As an example of the unique decay topologies generated by these models, we consider pair produc- 
tion of new exotic heavy quarks, pp QQ. The leading production mechanism is via QCD 

pp^ sQ* sQ + X ^qXsQ + X ^qxqllx + X (48) 

where X represents the beam remnant and other possible hadronic activity. The first sQ decays 
via sQ* — )• g X and the second decay to sQ qllx which is a primary decay chain of study. The 
charged leptons are I = e, fj,. The signal is for two isolated leptons, two light quark jets and large 
amounts of We take a mass spectrum of 

mq = 700 GeV = 650 GeV m^L = 300 GeV = 100 GeV my = 100 GeV (49) 
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jet 

~sQ I 7l dm 

Figure 6: The topology of the primary decay chain. Here DM is the x particle. 

The topology of the primary decay chain is shown in Figure [H The partial decay widths for the 
sQ is 




_ ^ ^ A|Mq I Ml 

M2 ^ levr V M2 



(50) 

'Q V ^"Q 

In the analysis, for simplicity, we set all of the Yukawa couplings to Ai = A2 = A. We assume a 
100% branching fraction of L — )• / sL and sL — )• / x- this model the gauge boson decays at one 
loop. sL is the lightest partner; thus, fastest the decay rate goes as 

^ 167r2 M8M4 ' 

A. 

where m, and M are the masses of the gauge boson, dark matter and sL, respectively. Here 
A is the coupling between the sL, x ai^d the SM lepton. We take A to have a technically natural 
value of A ~ 0.001 so the gauge boson is long-lived. With the masses given in equation HOj we 
have a lifetime of about 10~^ seconds. It should be noted that long-lived particles take about 
~ 0(1) X 10~^ seconds to transverse the larger ATLAS detector. Thus, these gauge bosons will 
register as missing energy. Even though the coupling is so small, the decay chain proceeds because 
of the branching fractions. Finally, the signal is generated with the new gauge bosons, V, being 
emitted from the decay chain in Figure [H We list the topologies generated to order a in Figures [82l 
[86] in Appendix B. In sections [2] and [Sj we have discussed the invariant mass distributions for dark 
matter stabilized with a Z2 or Z3 stabilization symmetry with the virtual particles, respectively, 
off- or on-shell. The present model presents a "hybrid" between the two pictures. This is because 
emitting the long-lived gauge boson forces part of the decay chain off-shell. Emitting the new gauge 
boson also causes these diagram to be suppressed because the virtual particles in the decay chain 
must go off-shell. Because there are three new gauge bosons, the overall off-shell suppression is 
enhanced by a multiplicity factor for each boson emitted. 

4.1.2 Extracting the Signal 

To get an estimate on the signal we first impose basic acceptance cuts which are consistent with 
ATLAS and CMS studies of on- as well as off-shell SUSY decay chains. [131 E] We require 

\r]i\<2.5, |??j|<2.5, (52) 

ARii > 0.3, ARij, ARjj > 0.4. (53) 
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where I and j are lepton and jets, rja is the pseudorapidity of particle a. ARab is defined as 
ARab = \J drj^^^ + d(j)'^i^ where drjab and dcl)ab is the difference in the pseudorapidity and transverse 
angle of the detector between particles "a" and "6." As described in the previous section, our 
example signal has two leptons and jets as the SM final states. The primary SM background for 
this signal is It decays into two leptons. Additional backgrounds include QCD, W+ jets and Z+ 
jets events. ATLAS and CMS places additional cuts to reduce this as well as other SM backgrounds. 
The model allows same-sign or opposite-sign dileptons in the final state. Since the purpose of this 
section is to see the effect of the detector cuts on our signal, we choose the more conservative 
opposite-sign dilepton cuts. We adopt cuts consistent with both collaborations by requiring 

1. Two leptons with pT > 20 GeV 

2. At least one leading jet with px > 100 GeV and subleading jets with px > 50 GeV 

3. It > 100 GeV and ^ > 0.2 M^s 

4. Transverse sphericity St > 0.2. 
Here the missing energy (fjr) is defined as 

$T=iT = - ^PiT (54) 

i 

and i runs over the transverse momentum, pT, of the visible final state particles in the event. The 
effective transverse mass, Meg, is defined as 

M,Q = Y,EiT + f}T (55) 

i 

where the sum runs over the measured transverse energy, Et, from the visible particles in the event. 
Finally the transverse sphericity {St) is defined as 

St = (56) 
where Ai^2 are the eigenvalues of the 2x2 sphericity tensor 

where k runs over the number of final state jets and leptons. The other indices, i and j, run over the 
Pt components of each particle. Sij is evaluated for the final states with rj < 2.5 and pT > 20 GeV. 
St ^ 1 for approximately spherical events; QCD events are usually back-to-back with St ~ 0. 
Generally, because our signal has three dark matter candidates per event, these cuts could be 
optimized with larger ^ cuts. For direct comparison with ATLAS and CMS, we simulated our 
signal with the cuts above. The ATLAS collaboration [13] finds the following backgrounds for 1 
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Background 


Events (1 fb-^) 


ti 


81.5 


W+ jets 


1.97 


Z+ jets 


1.20 


QCD 





Total SM 


84.67 



Table 2: SM backgrounds as computed by [13]. 



fb~^ of integrated luminosity (see Table 2). In addition to these backgrounds, we have an additional 
irreducible background when the Z2-like signal process, equation SHI emits Z boson which decay 
invisibly. The invisible branching for Z bosons into neutrinos is 20%. |25j Finally, in our analysis we 
simulate calorimetry responses for the energy measurements by adopting Gaussian smearing [13] 
with the following parameters. 
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4.1.3 Results 
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Figure 7: Dilepton invariant mass (left panel) and invariant mass squared (right panel) distributions 
for the topology in Figure [H The cuts described in section [1.1.21 are applied. 



We ran our Monte Carlo for the LHC at 14 TeV center-of-mass energy for the signal process 
in Figure m We used CTEQ 4M parton distribution functions |26j ; and, all results are presented 
at parton level. Us is computed at two-loop level. At zero order in the SU{2)£) gauge coupling, 
the model admits "Z2-like" topologies. We show the kinematic edge resulting from this topology 
in Figure [71 We also include SM the dominant tt as well as the irreducible Z ^ vv backgrounds. 
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Figure 8: Dilepton invariant mass (left panel) and invariant mass squared (right panel) distributions 
but with the first order corrections from the long-lived gauge bosons. The left panel features the 
double kinematic edge. The edges are roughly separated by the 100 GeV gauge boson mass. The 
Z2-like signal is also plotted for comparison as well as the backgrounds described in the figure 
above. All of the topologies generated by the long-lived gauge bosons are listed and individually 
plotted in Appendix B. 



To order a in the SU{2)£) coupling, we have six additional diagrams which generate corrections. 
For completeness, in appendix B, we list all of the different topologies and plot each correction 
before interfering the diagram to get Figure [HI Each diagram listed in the appendix is instructive 
for the shape and position for each kinematic edge. The kinematic cuts listed in section [4.1.21 are 
taken; the shapes of the distribution are generally preserved under the cuts. The total irreducible 
background events from Z ^ vv and the dileption tt channel for 100 fb~^ is 

Bz^vv = 98.5 Bit = 56630 (59) 

Additionally total signal events (Figure [8|) for 100 fb^^ 

S = 4440 (60) 

which generates the following signal-to-background ratio and statistical significance for signal ob- 
servability 

SjB = 0.08 = 18.6 (61) 

4.2 Warped GUT 

We present another very well-motivated Z3 model: for more details, see the original references 
[17^ I18j . This model is based on the framework of a warped extra dimension with the SM fields 
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propagating in it which can address both the Planck-weak and flavor hierarchy problems of the SM 
|27j . In a grand unified theory (GUT) model within this framework, it was shown that 

• a viable DM particle can emerge a spin-off oi suppressing proton decay. 
Moreover, 

• Z3 (rather than Z2) as the symmetry stabilizing DM arises naturally (due to combination of 
SM color and baryon number quantum numbers). 

It turns out that the SM particles resulting from the decay of mother particles in this model are 
mostly top quarks and Vl^'s. We defer an analysis of the reconstruction of these decay chains to 
future work. Here we simply give a summary of this model and the relevant LHC signals. 

4.2.1 Basic framew^ork 

In this framework, the SM particles are identified as zero-modes of 5D fields, whereas the heavier 
modes (i.e., non-trivial excitations of the SM particles in the extra dimension) are denoted by 
Kaluza-Klein or KK particles and constitute the new physics. A few TeV KK mass scale can be 
consistent with electroweak and flavor precision tests, at the same time avoiding at least a severe 
fine-tuning of the weak scale |27j . 

An extension of the SM gauge group in the bulk to GUT gauge group is motivated by precision 
unification of gauge couplings and explanation of quantization of hypercharge. In more detail, the 
extra/non-SM 5D gauge bosons (denoted generically by X) do not have zero-modes by, for example, 
an appropriate choice of boundary conditions. However, these gauge bosons still have KK modes 
with a few TeFmass (i.e., same as SM gauge KK modes), instead of usual mass of ~ O (lO^^) GeV 
in 4Z)-like GUTs. The fermions follow a different story as follows. The 5D fermion fields must of 
course form complete GUT multiplets. Usually, an entire SM generation fits in such a complete 
multiplet(s), for example, 16 for 50(10) GUT group, i.e., quark-lepton unification is incorporated. 
However, if we attempt to identify SM fermions of one generation as zero-modes of a complete 5D 
GUT multiplet, then it turns out that we will get too fast proton decay via exchange of X between 
SM quarks and leptons - again, with a few TeV massif 

4.2.2 Split fermion multiplets 

Fortunately, the breaking of GUT gauge group down to SM gauge group by boundary condition 
allows "split" fermion multiplets (just like for gauge bosons) as follows. We can choose bound- 
ary conditions such that one 5D multiplet (labeled "quark" multiplet) has zero-mode only in its 

^■^It turns out that the couphngs of X to SM quarks and leptons are suppressed - roughly by powers of SM Yukawa 
couplings - due to the nature of the profiles in the extra dimension of the various particles, but this effect is not 
sufficient to allow a few TeV mass for X to be consistent with proton decay. 
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quark component, with the lepton-hke component having only KK mode and vice versa for an- 
other multiplet (labeled "lepton" multiplet'). Thus SM quarks and leptons originate from different 
5D multiplets, avoiding exchange of X gauge bosons between SM quarks and leptons since such 
exchange can only couple fermions (whether zero or KK-modes) within the same 5D fermion multi- 
plet. In spite of this "loss" of quark-lepton unification, the explanation of hypercharge quantization 
is still maintained since SM quark must still be part of a complete GUT multiplet. In fact, such 
splitting of fermion multiplets results in precision unification of couplings in the model where the 
GUT group is broken down to the SM on the Planck brane [25]. The reason for the modification 
relative to the running in the SM (and hence the improvement in the unification) is the different 
profiles for quarks and leptons, especially within the third generation. 

4.2.3 Additional symmetry for proton stability 

It turns out that to maintain this suppression of proton decay at higher orders, we have to impose 
an extra symmetry, for example, a gauged U{1)b [commuting with the GUT group] in the bulk as 
follows. The entire quark multiplet is assigned B = 1/3 (i.e., that of the zero- mode contained in 
this multiplet). Thus lepton-like states from this multiplet are "exotic" in the sense that they have 
B = 1/3. Similarly, the entire lepton multiplet is assigned B = 0, giving exotic quark- like states 
(i.e., with B = 0). X^s are also exotic since they are colored, but have B = 0. The exoticness 
is especially striking since these states cannot decay into purely SM: explicitly, they are charged 
under a symmetry 

$ ^ e2-(^-^)cI> (62) 

(where a, a are the number of color, anti-color indices on a field <I>), under which SM is neu- 
tralH Thus the lightest Za-charged particle (dubbed "LZP") is stable (others Zs-charged particles 
produced in colliders or in the early, hot universe decay into it) and a potential DM candidate, 
depending on its couplings. 

4.2.4 Who's the LZP/DM? 

In the model with GUT group broken to the SM on the Planck brane which was the focus in 
references [IZIIIH], it turns out that, due to profile of tji (in turn, based on heaviness of top quark 
and constraint from shift in Zbb coupling), its exotic GUT partners are lighter that typical KK 
scale (say, mass of gauge KK modes). Thus, it is very likely that LZP resides in this multiplet. In 

^*In more detail, U{1)b has to be broken to avoid zero-mode gauge boson. We break it on the Planck brane so 
that 4-fermion operators giving proton decay [i.e., violating U{1)b] can only arise on the Planck brane where they 
are adequately suppressed. However, U{1)b still cannot be broken arbitrarily, i.e., a subgroup of U{1)b must still 
be preserved in order to forbid (mass) mixing of lepton and quark multiplets on the Planck brane which will lead to 
rapid proton decay - for example, we require that the scalar vev which breaks U{1)b has B = integer in which case 
the above Zs symmetry is still preserved, even if U{1)b is broken. 
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particular, if the GUT group is 50(10), then there is a GUT partner of tji with quantum numbers 
of right-handed (RH) neutrino, but with B = 1/3 (denoted by u'j^ and its LH Dirac partner, denoted 
by v'j^). It has been shown that if this z^' is the LZF0, then it is a good DM candidate, in the sense 
that, in some regions of parameter space, it has the correct rehc density upon thermal freeze-out 
in the early universe: see Fig. 5 of reference |18] and Figs. 3 and 4 of reference |29] (the latter 
reference studies a modified version of the model outlined here) and related discussion. Similarly, 
the constraints from direct detection of DM can be satisfied: see Fig. 7 of reference [18] and Fig. 
3 of reference and related discussion. Other GUT partners of tR are then heavier than i/' , but 
they can still be lighter than SM gauge KK modes. And, GUT partners of other SM particles, 
including X-type gauge bosons, are as heavy as SM gauge KK modes. 

4.2.5 DM partner at the LHC 

As usual, in order to test this idea at colliders, we consider producing the (other than LZP) Z^- 
charged particles at colliders (of course, these must be produce in pairs) and observe their decays 
into SM particles and LZP. Since colored and lightest such particles will have largest cross-section at 
the LHC, a good candidate for such a study is the GUT partner of tji with {t,b)L quantum numbers, 
denoted by {t'j^, h'^) [and it's conjugate by it'^, ^^'j]!^ The two states t' and h' are degenerate before 



EWSB, but will be split afterward. 

We focus here on b' due to the interesting features in its decay channels as shown in figure 
[91 Xs, the SU{2)l doublet X, Y and another SU{2)l doublet X', Y' are beyond SM gauge 
bosons of 50(10), with electric charges 2/3, 4/3 and 1/3, and 2/3 and 1/3, respectively0 First 
of all, the 1st decay chain of b' into two DM {v'^ and tW does have the topology needed to 
give a cusp in the invariant mass of SM/visible particles (of course assuming on-shell intermediate 
particles; see comment below). Second, the 2nd process above involves the same final state, but 
with a different topology and thus is relevant for obtaining multiple edges (again, with on-shell 
intermediate particles), even with two DM in final state. However, it is clear that the intermediate 
particles (X-type gauge bosons) in the above processes are actually OjQ^-shell so that these features 
are not so useful for us here. Finally, while there does not seem to be a decay of b' into tW and one 
DM which would be relevant for obtaining double edge with off-shell intermediate particles (when 
combined with the 1st and 2nd processes above with two DM), there is a decay to tWW and one 
DM (i.e., with an "extra" W) as in the 3rd process above which might play this role. Thus, decay 
of b' will exhibit a double edge due to presence of one DM and two DM in final state [along the 
lines of the discussion of Eq. (j2.1.2p ]: M^/ — 2 M^i — mt — mw (for 1st and 2nd processes) and 
Miji — Myi — 2 rrit — mw (for 3rd process). 



^^At leading order, all the GUT partners of are degenerate, but higher-order efTects can split them. 
^''Recall that the {t,b)L and conjugate of tn are contained in the same representation, namely, 16 of 5*0(10) so 
that the {t'i^,h'i) states with B — —1/3 are indeed exotic. 

^'^Note that only the 1st decay channel is shown in Fig. 10 of of reference [18] . 
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Figure 9: Different possible decay cliains for b' in the scenario of [18j. Note the appearance of two 
DM states in two of the possible decay chains, whereas only one DM particle comes out of the third 
possibility. 



Note that, in a non-minimal model, for example if SO(IO) is broken on the TeV brane instead 
of the Planck brane, the X's bosons might be lighter so that the intermediate particles in decay 
chains similar to those above might be on-shell. However, then we might as well study production 
of the lighter X's (instead of the exotic fermions) which are also colored and which will decay into 
LZP via ojff-shell GUT partners of SM fermions. 



5 Conclusions and Outlook 

Many extensions of the SM contain stable WIMPs which can be viable DM candidates. Most 
of these models stabilize the DM with a parity {Z2) symmetry. However, in the spirit that all 
fundamental particles in nature are defined by how they transform under different symmetries, 
these models actually correspond to one type of candidate! On the other hand, any continuous or 
discrete global symmetry can be adopted for DM stabilization; and, DM candidates stabilized by 
a parity symmetry and, e.g., a Z3 symmetry are different. 

This possibility of other than Z2 symmetries stabilizing DM is more than just of academic 
interest; the nature of the stabilization symmetry has important implications on collider searches 
for DM. At colliders other particles (heavier than DM) which are charged under the same symmetry 
which stabilizes the DM candidate(s) can be produced, decaying into DM and SM particles. Such 
events will generate decay topologies and modes that are determined, in part, by the stabilization 
symmetry. Thus the analysis of such "missing energy" signals will present a hint not only for the 
existence of the DM but also the nature of its stabilization symmetry. 
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For example, the decay of a single such mother particle can contain one or two DM in the case 
of a symmetry, but only one in the case of Z2 symmetry. We showed that in many cases simple 
kinematic observables, such as invariant mass distributions of the visible particles of such decay 
chain, could then characterize the stabilizing symmetry. Specifically, when a mother particle decay 
via off-shell intermediate states into the same visible particles along with one and two DM (for 
the case of symmetry), it may be possible to observe a double-edge in the distribution of these 
visible particles (vs. single edge for Z2 symmetry). In fact, the difference between the location of 
the edges will be a direct measure of the mass of the dark matter particle for Z^ models. On the 
other hand, when the intermediate particles are on-shell, we also pointed out the possibility of a 
very distinctive feature appearing in the invariant mass distribution of two visible particles in the 
case of Z3 symmetry: a cusp dividing the distribution into two regions. This happens when two 
DM particles emerge from the same chain, with one of these DM particles being situated in-between 
the two SM particles. 

Signal Fakes: We point out that models with parity stabilized dark matter can naively fake 
double kinematic edges and cusps. An example of a faked cusp comes from a decay chain with 
on-shell intermediate particles and involving a SM neutrino(s) in the final state. If a neutrino is 
emitted between the two other, visible SM particles, then we obtain the new topology considered 
above so that the invariant mass distribution of the two visible SM particles can give a cusp. As 
another example, a double kinematic edge also comes from two different mother particles charged 
under Z2 which could be produced in separate events [30], but with same visible decay productJ^ 
Also, in this paper, we considered decays of mother particle into DM and SM particles occurring 
inside the detector. However, a long-lived, parity odd particle in a Z2 model that decays outside 
of the detector can fake our Z3 signal. For example, a single mother can decay into such a particle 
and SM particles. The same mother can also decay to the same SM final state and the DM. A 
combination of these two decay chains can generate a double edge. A closely related scenario is 
that there are actually two (absolutely) stable DM particles in a Z2 model [3l], again giving double 
edges from two decay chains (along the lines just mentioned). Finally, if there are more than two 
SM particles involved in the decay of a single mother, the analysis of the cusp in the invariant mass 
distribution will be more complicated, but nevertheless we expect that in general the same type of 
arguments made here should be able to be implemented in this case. As an example, a possible 
pseudo-faking situation will arise when three visible particles are emitted in a decay cascade in the 
Z2 case. However, we mentioned how this type of faking can be resolved by considering all possible 
pair invariant masses. 

Future Considerations: Of course, in any given event, there will be two such mother par- 
ticles presenile (three mothers is a possibility in Z3 case, but it is phase-space suppressed). The 

^^Note that, assuming pair production, multiple mothers in Z2 cannot fake different edges in two distributions that 
we discussed in section [521 

^^However, studying the decay of only one mother (as we have mostly done in this paper) is still relevant, for 
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assumption of a Zi symmetry thus points to an eventual emergence of two invisible particles for 
every new physics event in the collider. On the other hand, models where dark matter is stabilized 
with, e.g. Z3 symmetry, can have two, three, or four dark matter candidates in an event. In an 
ongoing work, we construct a variant of the and rriTX variables [3] to use the information on 
both decay chains in the full event. The goals of such techniques are to establish that the DM 
is stabilized by a Z3 symmetry in other situations which were not discussed here, e.g., when a 
cusp-like topology is absent, and to better eliminate the fakes of Z3 signal by Z2 models. We also 
are studying other techniques to eliminate the fakes described above. 

In all, we emphasize that parity symmetries are not necessarily the only way to stabilize the 
DM and we showed (via a few example cases) that models with other stabilization symmetries 
generally have testable consequences at the LHC, i.e., can potentially be distinguished from the 
parity case. The reader should regard this work as a first step into a more complete study of beyond 
Z2 stabilized dark matter scenarios. 

Acknowledgments 

We would like to thank I. Hinchliffe, A. Katz, G. Servant, M.D. Shapiro, R. Sundrum and B. Tweedie 
for valuable discussions. K.A. was supported in part by NSF grant No. PHY-0652363. D.W. is 
supported in part by a University of California Presidential Fellowship. 

References 

[1] See, for example, G. Bertone, D. Hooper and J. Silk, Phys. Kept. 405, 279 (2005) 
| arXiv:hep -ph/ 0404175l . 

[2] G. Hinshaw et al. [WMAP Collaboration], |aj2£iYi0803,0732J [astro-ph]. 

[3] E. W. Kolb and M. S. Turner, Front. Phys. 69, 1 (1990). 

[4] C. G. Lester and D. J. Summers, Phys. Lett. B 463, 99 (1999) [arXiv:hep-p h79906349 1 ; 
A. J. Barr, C. G. Lester, M. A. Parker, B. C. Ahanach and P. Richardson, JHEP 0303, 
045 (2003) [arXiv:hep-ph/02 08214|; A. Barr, C. Lester and P. Stephens, J. Phys. G 29, 2343 
(2003) |arXiv:hep -ph/0 304226| ; W. S. Cho, K. Choi, Y. G. Kim and C. B. Park, Phys. Rev. 
Lett. 100, 171801 (2008) [arXiv:0709.0288l [hep-ph]]; B. Gripaios, JHEP 0802, 053 (2008) 
|arXiv:0709.2740] [hep-oh]]: A. J. Barr, B. Gripaios and C. G. Lester, JHEP 0802, 014 (2008) 
|arXiv:0711.4 008 [hep-ph]]; W. S. Cho, K. Choi, Y. G. Kim and C. B. Park, JHEP 0802, 035 
(2008) [arXiv:0711.4526l [hep-ph]]. 

example, if this decay involves a "clean" (SM) final state (vs. a dirtier one from the other mother). Of course, we 
also need to determine experimentally which SM particles in the event came from this decay chain, which is possible, 
for example, if this mother is boosted in the laboratory frame so that the decay products are roughly coUimated. 



32 



[5] H. C. Cheng, D. Engelhardt, J. F. Gunion, Z. Han and B. McElrath, Phys. Rev. Lett. 100, 
252001 (2008) |arXiv:0802.4290l [heo-ph]] : H. C. Cheng, J. F. Gunion, Z. Han and B. McElrath, 
Phys. Rev. D 80, 035020 (2009) |arXiv:0905. 13441 [hep-ph]]. 

[6] K. T. Matchev, F. Moortgat, L. Pape and M. Park, JHEP 0908, 104 (2009) tarXiv:0906.2417] 
[hep-ph]]; P. Konar, K. Kong, K. T. Matchev and M. Park. larXiv:0911.4126l [hep-ph]. 

[7] T. Han, I. W. Kim and J. Song. " arXiv:0906.5009l [hep-ph]. 

[8] G. Jungman, M. Kamionkowski and K. Griest, Phys. Rept. 267, 195 (1996) 
|arXiv:hep- ph/9506380] . 

[9] H. S. Lee, Phys. Lett. B 663, 255 (2008) [ar Xiv:0802.050 6 [hep-ph]]. 

[10] H. C. Cheng and L Low, JHEP 0309, 051 (2003) and JHEP 0408, 061 (2004) 
|arXiv:hep-ph/0405243| . 

[11] G. Servant and T. M. P. Tait, Nucl. Phys. B 650, 391 (2003); H. C. Cheng, J. L. Feng and 
K. T. Matchev, Phys. Rev. Lett. 89, 211301 (2002) [ arXiv:hep-ph/0207125] . 

[12] K. Agashe, A. Falkowski, L Low and G. Servant, JHEP 0804, 027 (2008). 

[13] The ATLAS Cohaboration, CERN-LHCC-99-015. 

[14] The CMS Cohaboration, J. Phys. G: Nucl. Part. Phys. 34 995-1579 (2007). 

[15] D. G. E. Walker. [a?Xiv:0907.3146l [hep-ph]. 

[16] D. G. E. Walker. [a?Xiv:0907.3142l [hep-ph]. 

[17] K. Agashe and G. Servant, Phys. Rev. Lett. 93, 231805 (2004) [arXiv:hep-ph/0403143| . 

[18] K. Agashe and G. Servant, JCAP 0502, 002 (2005) ^arXiv:hep-ph/0411254| . 

[19] E. Ma, Phys. Lett. B 662, 49 (2008) [arXiv:0708.3371 [hep-ph]]. 

[20] E. Byckling and K. Kajantie, Particle Kinematics (John Wiley & Sons, 1973). 

[21] D. J. Miller, P. Osland and A. R. Raklev, JHEP 0603, 034 (2006) [arXi v:hep-ph/051 0356] . 

[22] S. Kraml and A. R. Raklev, Phys. Rev. D 73, 075002 (2006) [arXiv:hep-ph/0512284| . 

[23] L. T. Wang and I. Yavin, JHEP 0704, 032 (2007) [arXiv:hep-ph/0605296| . 

[24] J. AlwaU et al, JHEP 0709, 028 (2007) (arXiv:0706.2334l [hep-ph]]. 

[25] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667, 1 (2008). 

33 



[26] H. L. Lai et al, Phys. Rev. D 55, 1280 (1997) |arXiv:hep-ph/9606399]. 

[27] For a review and the original references, see H. Davoudiasl, S. Gopalakrishna, E. Ponton and 
J. Santiago, arXiv:0908.1968 [hep-ph]. 

[28] K. Agashe, R. Contino and R. Sundrum, Phys. Rev. Lett. 95, 171804 (2005) 
|arXiv:hep-p h /0502222'i . 

[29] K. Agashe, K. Blum, S. J. Lee and G. Perez, Phys. Rev. D 81, 075012 (2010) [arXiv:0912.3070l 
[hep-ph]]. 

[30] For a supersymmetric example see, H. Baer, A. Mustafayev, E. K. Park and X. Tata, JHEP 
0805, 058 (2008) [ar Xiv:0802.3384. [heo-ph]]. 

[31] For some recent studies, see, for example, K. M. Zurek, Phys. Rev. D 79, 115002 (2009) 
[arXiv:08n.4 429 [hep-ph]]; F. Chen, J. M. Cline and A. R. Frey, Phys. Rev. D 80, 083516 
(2009) [arXiv:0907.4746t [hep-ph]]: D. Feldman, Z. Liu, P. Nath and G. Peim, Phys. Rev. D 81, 
095017 (2010) [arXiv: 1004 .0649 [hep-ph]]; M. Cirelh and J. M. Cline. larXiv:1005.1779l [hep-ph]. 

Appendix 

A. The distribution for the new topology 

Most of the intermediate steps in the derivation of the cusp in Eq. (j28p are similar to the analysis 
in reference |21j of the reaction in Eq. (|2ip . except for the fact that a DM (i.e., massive) particle 
is situated in-between two SM particles in the new topology (See Eq. (|24p ). Based on the algebra 
and the notations found in reference [21], we will derive a few useful relations. 

Basically, the invariant mass formed by the two SM particles in this topology is given by 

= {Pc + Pa? = 2E^Ea{l - COS ^ea) (63) 

where 9ca is the opening angle between two visible particles. Note that this relation is always valid 
in any frame so that we can rewrite the above relation as 

mL = 2i?(^)4^)(l-cosC^). (64) 

Here and henceforth the (particle) superscripts on 0's (in this case B) imply that those angles are 
measured in the rest frame of the corresponding particle. Using energy-momentum conservation, 
we can easily obtain the energies for a, DM, and c, which are measured in the rest frame of particle 
B. 
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^ !^ (65) 

<j - "^^"7^;""- (66) 

^(B) ^ {ml - ml)ms ^^^^ 

m%+ml;- - AV2(m2,, m|, ml^) cos 6*^ 

Inserting these relations into Eq. (j64p . we obtain 

"^ca = 2 ^^-^ Tg^ ^ ^(l-cos^(f)). (68) 

We easily see that the maximum of m^^ occurs when cos 0^ j-jj^ = 1 and cos Oca = — 1- We want to 
express the invariant mass ruca in terms of variables which have flat distributions: this is the case 
for cos9ca \ but not for cos^^^^jyj. So, we need to express cos^^^^j^ in terms of cos^^^^j^ (i.e., the 
same angle in the rest frame of particle C) for which the distribution is also flat. This relation can 
be found by calculating "t-^dm frames of particle C and B: 



2 _ 2 

"^c DM — ""^DM 



+ 2EPe^^^ - 2EP J{E^^^Y - ml^ cos of^u (69) 



= ml^ + 2Ei^)E^^^-2Ei^)^ {E^^^f -ml^cosei^'^M (70) 
Again, the energy-momentum conservation in the rest frame of C gives the following relations: 

- ^ ^ 

eP = H!^^ (72) 
2mc ^ ' 



Substitution of Edm and Ec in the rest frame of C and B into Eq. ()69p and Eq. ()70p gives the 
relation between cos^jj^]^ and cos^^^j^: 



DM 



m2. + m| -m2)j^ - AV2(TO2,,rn|,m2)j^)cos0^DM \ 
Next, we introduce the variables u and v: 

^_l^co^^^ ^_ 1-cosg^f^ ^^^^ 

and using Eq. ([73j) . we express m^^ in terms of u and v: 

ml, = {mZn\'^-o^u)v (75) 



(73) 
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where 



[nT'ca ) 



2(m|, — m^(j){m\ — m.\) 



''DM 



(76) 



(77) 



Note that the differential distributions for u and v (jd < u,v < 1) are also flat: 

where 6{x) is the usual step function. Replacing u and u by u and m^^ by using Eq. (f75|) gives the 
differential distribution 

1 d^T 



r dudmi^ 



e{u) 



(m™a^)2(l - au) J (m™a^)2(l - on) 



(78) 



where a "top-hat" function 9[x) = 9{x)9{l — x). The next step is to integrate over u to find the 



distribution in m^^: 



1 d^T 



1 a^r 



du 



du 



-du 



(m--)2(l - au) 



for < ruca < ^m^^) where 



tin 



Max 1, 



1 



m;: 



(^max)2 



(79) 



(80) 



Now the above integral is easy to evaluate, and we finally obtain the distribution which was given 
earlier in Eq. ([28]) : 



1 d^T 



1 



^In^ 



for 



< < \/r 



am;: 



(81) 



for Vl — amfa"^ < rrica < 



m" 



B. Signal Topologies from Section 14.1.11 

In addition to the decay chain in figure [6] for the model presented in section [5?T1 there are additional 
corrections by the long-lived SU (2) d gauge bosons. The masses are listed in Eq. |39l As described 
above, the virtual particles in the decay chain go slightly off-shell when emitting the new gauge 
boson; however, because there are three of them, the additional multiplicity factor helps to amelio- 
rate this suppression. In this appendix for each topology, we plot the invariant mass distributions 
with the cuts in section 14.1.21 (see Figure [TOlllSp . 
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Figures [82]|85} Order a corrections by the SU (2) £> gauge bosons to the decay chain in Figure [H 
See the model in section 14.11 
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Figure (86) Order a corrections by the SU (2) ^ gauge bosons to the decay chain in Figure [6l See 
the model in section 14.11 
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Figure 11: The plots corresponding to the topology of Figure 
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Figure 12: The plots corresponding to the topology of Figure 
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Figure 13: The plots corresponding to the topology of Figure 
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Figure 15: The plots corresponding to the topology of Figure 
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